Calculations of Magnetic Exchange Interactions in Mott— Hubbard Systems 
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An efficient method to compute magnetic exchange interactions in systems with strong correlations 
is introduced. It is based on a magnetic force theorem which evaluates linear response due to 
rotations of magnetic moments and uses a novel spectral density functional framework combining our 
exact diagonalization based dynamical mean field and local density functional theories. Applications 
to spin waves and magnetic transition temperatures of 3d metal mono-oxides as well as high-T c 
superconductors are in good agreement with experiment. 
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Obtaining a quantitative theory of magnetic materials 
spanning from itinerant to atomic limit, above and be- 
low their temperatures of magnetic ordering has been a 
theoretical challenge for many years 0. By now itiner- 
ant magnets are well described by local spin density ap- 
proximation (L5DA) of density functional theory (DFT) 
2], where methodologies based on spin-spiral frozen- 
magnon technique , the use of magnetic force theorem 
via an evaluation of linear response due to rotations 
of magnetic moments as well as spin dynamics calcula- 
tions [5j have allowed to access a great variety of physical 
properties such as spin wave spectra , magnetic ordering 
temperatures, interatomic exchange constants, dynami- 
cal susceptibilities, etc [ft m. M. M. llC| . 

However, there is large class of systems where cal- 
culations of exchange interactions is still a challeng- 
ing theoretical problem. These are strongly correlated 
systems like high-T c superconductors or atomic mag- 
nets where the on-site Coulomb interaction U is com- 
parable or larger than the bandwidth. In cases where 
magnetic ordering and/or lattice distortions result in 
a non-de gen erate equilibrium state, techniques such as 
LDA+U ITU or GW0 have been applied to describe 
spectroscopy, magnetic moments, and even spin wave 
spectra of systems such as MnO However, in gen- 

eral, excitation spectra of strongly correlated systems are 
not representable by single Slater determinants and show 
such features as atomic multiplets 0], Zhang-Rice sin- 
glets jl5| . Kondo resonances, etc. In Mott-Hubbard in- 
sulators the energy gap is much larger than the magnetic 
ordering temperature above which a local moment regime 
takes place, i.e. the system becomes paramagnetic but 
remains insulating. These properties cannot be accessed 
either by static mean field approaches such as LS'DA or 
LDA+U or by perturbation theory over the Coulomb in- 
teraction such as GW. While versions of static as 
well as dynamic 01 coherent potential approximations 
have been introduced in the past to access disordered lo- 
cal moment regime, developing a generalized framework 
having a capability to compute exchange interactions in 
both itinerant and atomic limits as well as in many in- 



termediated cases would open new opportunities in com- 
putational design of new magnetic materials. 

In the present work we develop a novel approach which 
is capable to deal with this problem. Our method is 
based on LDA+Dynamical Mean Field Theory (DMFT) 
[l8| . a recently developed electronic structure method 
describing correctly both itinerant and atomic limits 
and accessing ordered and disordered moment regimes 
on equal footing. This is achieved by treating corre- 
lated electrons with frequency dependent self-energies 
deduced from solving corresponding Anderson impurity 
problem (AIM) subjected to a self-consistency condition. 
Conveniently formulated using a spectral density func- 
tional 0, this LDA+DMFT method incorporates re- 
alistic band structure effects and has already helped to 
solve several long-standing problems, e.g. temperature- 
dependent magnetism of Fe and Ni[l9j, volume collapse 
in Cej^fji and Puji^], as well as electronic structure of 
doped Mott insulators^. 

In order to deduce exchange constants for general wave 
vector q we utilize a linear response based magnetic 
force theorem 0,|i3]and demonstrate the accuracy of our 
method on several Mott-Hubbard insulators including 
late transition metal oxides as well as parent materials of 
high-T c superconductors. To solve the impurity problem 
which is the nucleus of the dynamical mean field algo- 
rithm we apply a newly implemented cluster exact diag- 
onalization method and calculate self-consistently local 
Green functions, self-energies and static linear response 
functions. Both the deduced spin wave spectra as well as 
magnetic ordering temperatures evaluated using a Monte 
Carlo simulations of the mapped Heisenberg Hamiltoni- 
ans are found in good agreement with experiment. 

Our implementation is based on a most recent many- 
body band structure algorithm |24| which allows us to 
avoid computationally expensive solution of the Dyson 
equation [u> — Hq(\s.) — E(w)]G(k, u>) = 1 for the electronic 
Green function G(k,u>) at a large grid of frequencies ui. 
This is achieved by assuming a pole expansion for the 
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self-energy 



E(w) = E(oo) + Y/?^ - P ^ Vi 



(1) 



so that the entire problem is reduced to a matrix equation 
with an energy-independent Hamiltonian 



u - ff (k) - E(oo) V"+ 



V" 



P 



G(k,w) =/ 



(2) 



where I is the unit matrix and an auxiliary Green func- 
tion C/ a/ g(k, w) is a matrix in the space of poles, while the 
physical Green function G(k, ui) corresponds to the first 
element of 0(k, to). In Eq.Q weights , and poles Pi 
can be viewed as matrices which provide a best fit to real 



It is remarkable that the present procedure allows us to 
use an expression for the interatomic exchange constants 
similar to a linear response formula derived within DFT 
. Consider second-order change in the total energy re- 
lated to the rotations of the magnetic moments appeared 
at sites R+t and R' + t' of the lattice (here R are the lat- 
tice translations and r are the atoms in the basis). The 
local magnetic field B T at every atomic site r is approx- 
imately described by the values of the self-energy taken 
at to = oo [for example, B% = EjT (oo) — E^(oo)]. Thus, 
according to the magnetic force theorem which assumes 
a rigid rotation of atomic spin, a linear response theory 
expresses the interatomic exchange constants in the form 



tRt'R' 



EE 

q \tjj> 



JUj /k+qj' 
e kj ~ £k+qj' 



(lpkj\[(T X B T ] Q |'0k+qj')(V'k+qj'|[ cr X Br/^IV'kj) 



„iq(R-R') 



(3) 



Here, a is the Pauli matrix while the one-electron energy 
bands e^ and quasiparticle wave functions ipkj appear as 
the solutions of the equation J5J), when using quasiparti- 
cle representation for the Green function <?(k, u>) in the 
form. 



g(k,u) = Y 



u - e k j 



(4) 



While viewed non-interacting like, this formula indeed 
contains major information about many-body features in 
the excitation spectrum. In particular, multiplet transi- 
tions as well as delocalized parts of the electronic states 
are represented by separate "energy bands" ekj includ- 
ing its k dispersion which is borrowed from the non- 
interacting Hamiltonian iJo(k). Thus, genuine redistri- 
bution of spectral weight driven by the many-body in- 
teractions is correctly captured by the present method 
which will give an important feedback on the calculated 
exchange interactions. 

There are two essential approximations which are made 
to make the theory computationally tractable. As has 
been discussed recently [23], the magnetic force theorem 
can be introduced for a Lattinger-Ward functional which 
would involve calculations of full frequency dependent 
integrals between the self-energies and the Green func- 
tions. The present method utilizes (i) the Hartree Fock 
values for the local magnetic fields, and (ii) rational fit, 
i.e. Eq. 0), to the self-energy, which allows us to per- 
form all frequency sums analytically while retaining all 
major many-body multiplet features of the spectrum in 



the convenient linear response expression (j3J). 

To illustrate the method we consider several 
transition-metal oxides MnO, FeO, CoO, NiO as well as 
parent high-T c compound CaCu02- All these materials 
are antiferromagnetic insulators with an energy gap of 
a few eV and Neel temperatures Tjv of a few hundred 
K. Staying below Tjv it is well known that LSDA signif- 
icantly underestimates the band gap of MnO and NiO 
and fails to predict insulating character for FeO, CoO 
and CaCu02- The LDA+U corrects for these failures 
but needs to assume a symmetry breaking for FeO and 
CoO. It is clear that being a Hartree Fock approximation 
the LDA+U would converge to a single Slater determi- 
nant ground state, while in many cases either degener- 
acy of the latter or proximity of low-lying excited states 
needs to be included in statistical averagings for the one- 
electron Green functions. All static mean field theories 
would necessarily fail to describe paramagnetic insulat- 
ing behavior. On the other hand, the LDA+DMFT is a 
method valid both in ordered and local moment regimes. 
Here, we consider the d electrons of transition metal el- 
ements as strongly correlated thus requiring dynamical 
treatment using DMFT. The s and p electrons are as- 
sumed to be weakly correlated and well described by the 
LDA Hamiltonian Rlda (k) including the full potential 
terms of the linear muffin-tin orbital (LMTO) method 
[j^j ]. To consider relativistic effects, the spin-orbit cou- 
pling is taken into account in all cases. To obtain the 
one-particle potential Ho(k) — HldaQ*) — Vdc entering 
(|2*|) we subtract the double counting term Vdc as pre- 
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scribed by Ref. We use the experimental lattice 

structure for all materials. 

In order to solve the Anderson impurity problem we 
implement cluster exact diagonalization (ED) method. 
For transition metal mono-oxides the clusters are chosen 
to include d orbitals of transition metal ions hybridized 
with oxygen p orbitals in the octahedral environment. It 
has been known for many years that such cluster exact 
diagonalizations provide a good description of photoe- 
mission spectra in these materials |26t v^\ . For the 2D 
system, such as CaCu02 this is reduced to a Cu d or- 
bital surrounded by an oxygen square. This treatment 
allows us to capture both an effect of atomic multiplcts 
and of the Zhang-Rice singlet being the lowest lying 
excitation of undoped high-T c 's. During iterations to- 
wards self-consistency, the Anderson impurity problem 
is exact diagonalized each time with the positions of 3d 
levels, the d-p hybridization matrix elements as well as 
the O 2p levels are extracted from the HldaO^-)- The 
latter being a density functional is allowed to recompute 
and readjust the parameters of AIM. The values of the 
Coulomb interaction U as well as the Hund's exchange J 
were obtained earlier by the constrained LDA calculation 
|2sj | and kept fixed throughout the calculation. The AIM 
gives access to the frequency dependent self-energies for 
the d electrons which are then rationally approximated 
assuming three poles fit in the formula (JTJ. 

Our calculated ground state properties including mag- 
netic moments and energy gaps are found to be in good 
agreement with experiment both below and above mag- 
netic ordering temperatures. Below Tjv this result is in 
accord with previous LDA+U studies |2j| and our nu- 
merical results indeed show that dynamical correlations 
only marginally influence values of magnetic moments in 
cases such as MnO, NiO and CaCu02- On the other 
hand, for CoO and FeO, their t2 g bands are only par- 
tially occupied and the ground states become degener- 
ate. Therefore small value of spin-orbital coupling has 
a large effect on an appeared orbital moments which are 
evaluated to be 0.36 and 1.02 hb in our cluster exact di- 
agonalized LDA+DMFT calculations for FeO and CoO 
respectively. Note that the LDA+U would be capable 
recover the insulating character only by assuming sym- 
metry lowered orbitally ordered solution. 

We now discuss our predictions for the calculated Neel 
temperatures. Based on the self-consistently obtained 
local Green functions and self-energies we evaluate the 
interatomic exchange constants as the integral over the 
q space using (6,6,6) reciprocal lattice grid and uti- 
lize Monte-Carlo simulations Q of the correspondingly 
mapped Heisenberg Hamiltonians to find Tn. Our re- 
sults are given in Table I, where together with our most 
accurate exact diagonalization based simulations we also 
list the predictions using LSDA and less accurate im- 
purity solvers: single atom exact diagonalization known 
as Hubbard I solver 0], and the Hartree-Fock approx- 



imation to atomic self-energy known as LDA+U. It is 
well known that LSDA underestimates the energy gap, 
and as a result it significantly overestimates the mag- 
netic transition temperatures as seen in Table I. (Since for 
FeO, CoO and CaCu02, LSDA converges to a completely 
wrong metallic state, we omit quoting those predictions). 
The LDA+U method considers the effect of U in a static 
way, and this partially cures the shortcomings of LSDA 
making Tjv smaller but still much larger than the exper- 
imental ones. This conclusion is in agreement with the 
previously calculated results for MnO [13| . Hubbard I ap- 
proximation uses exact diagonalization for atomic d-shell 
and deduces frequency dependent self-energy, which fur- 
ther reduces Tjv for all transition-metal oxides. The best 
results are seen to be obtained by allowing d-electrons to 
fluctuate between the bath and impurity as prescribed by 
the cluster exact diagonalization calculation. This effect 
is missing in the Hubbard I solver and seriously affects 
the electronic structure of studied materials as, for ex- 
ample, for CaCu02 it redistributes the spectral weight 
by bringing such features as the well known Zhang-Rice 
singlet 15]. In an extreme situation, where the hybridiza- 
tion is much larger than the local Coulomb U, the mag- 
netization would eventually disappear due to strong fluc- 
tuation in the number of (i-electrons at the impurity site. 
So it is easy to understand why the cluster ED method 
gives smaller magnetic transition temperatures which are 
now closer to the experiment. 

It is interesting to discuss the physical reasons why Tjv 
decreases when going from NiO to MnO. These Mott- 
Hubbard insulators show almost atomic values of mag- 
netic moments M = 10 — n corresponding to d n configu- 
rations, which would under assumption of the same anti- 
ferromagnetic exchange constant Jaf mean that the or- 
dering temperatures should increase with increasing the 
moments. However, Jaf will decrease significantly due 
to the change in the lattice parameter. Also, during eval- 
uation of Tn we need to account for the quantum aver- 
aging for atomic spins directions which gives a prefactor 
S(S + 1)/S 2 deviating from 1 for small S. To sort out 
these effects we first performed a sample calculation for 
NiO with the expanded lattice constant of MnO. The 
Tjv for NiO has dropped from 519K to 327K in this case. 
Second, since Sxi = 1 while Smu = 5/2, the prefactor 
S(S + l)/S 2 would account for a 40% difference so that 
327x1. 4/2=229K is the Neel temperature that we need 
to compare with our predicted Tjv=172K for MnO. The 
residual discrepancy can be attributed to different ex- 
change splittings which also affects Jaf as pointed out 
earlier [3l| . 

We now discuss our calculated spin-wave dispersions 
along major symmetry directions in the Brillouin Zone. 
We illustrate this calculation on a case of NiO for which 
the magnon spectra have been measured long time ago 
[32I ]. Fig. 1 shows the results of our simulations using 
several levels of approximations. In accord with our pre- 
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TABLE I: Comparison between calculated using various ap- 
proximations and experimental magnetic transition tempera- 
tures (in K) in selected Mott-Hubbard systems. Hubbard I 
and Cluster ED denote the results of LDA+DMFT calcula- 
tions using Hubbard I and cluster exact diagonalization im- 
purity solver. 



LSDA LDA+U Hubbard I Cluster ED Exp. 



MnO 423 


240 


180 


172 


122 a 


FeO 


344 


297 


211 


198 a 


CoO 


407 


356 


300 


291 a 


NiO 965 


603 


542 


519 


523 a 


CaCu0 2 - 


765 


698 


602 


537 b 



a Ref. 29]; b Ref. 30]. 
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FIG. 1: Theoretical spin-wave dispersions for NiO calculated 
by LSDA, LDA+U, Hubbard I and Cluster Exact Diagonal- 
ization (ED) diagonalization impurity solvers in comparison 
with the experiment [32j. 

dictions for the Neel temperatures, the spin waves are se- 
riously overestimated by the LSDA theory but get closer 
to the experiment once correlations are taken into ac- 
count. The best accuracy is achieved when using the 
cluster exact diagonalization as the impurity solver in the 
LDA+DMFT calculation which demonstrates the impor- 
tance of many body redistribution of the spectral weight 
on the calculated exchange integrals. 

In summary, based on our newly implemented clus- 
ter exact diagonalization LDA+DMFT calculation, we 
presented a linear response method to calculate the ex- 
change interaction parameters of strongly correlated sys- 
tems valid as long as mapping of total energy functional 
to rigid spin based Heisenberg Hamiltonians makes sense. 
By using the rational interpolation for the self-energy, 
our approach is very efficient, and this has allowed us to 
describe quantitatively spin-wave dispersions and mag- 
netic transition temperatures of several realistic Mott- 
Hubbard insulators with many atoms per unit cell. Ap- 
plications to metallic systems are more challenging as 



they may need much larger clusters to account for such 
subtle effects as, e.g., the Kondo screening, and will be 
carried out in the future work. 
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